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Summary 


A theory  of  scattering  of  seismic  body  waves  by 
small  random  spacial  fluctuations  in  density  and  elastic 
parameters  in  an  otherwise  spherically  symmetrical  Earth 
model  is  developed.  It  is  shown  that  a primary  wave 
disturbance  of  either  P or  S type  travelling  through 
a slightly  irregular  solid  medium  will  generate  scattered 
waves  of  both  P and  S types.  Explicit  formulas  are 
derived  for  the  mean  square  amplitudes  of  P waves  scattered 
from  several  different  assumed  forms  of  primary  P wave. 

The  theory  assumes  that  the  primary  wave  may  be 
locally  approximated  by  a plane  wave  inside  each  part  of  the 
scattering  region  whose  size  is  comparable  with  the  mean 
size  of  the  irregularities  present  and  that  ordinary  ray 
theory  may  be  applied  to  calculate  the  travel  times  and 
amplitudes  of  both  primary  and  scattered  waves. 


The  theory  supports  the  writer's  earlier  suggestion 
that  observed  precursors  to  the  seismic  core  phase  PKIKP 
may  originate  by  scattering  from  the  phases  PKP^  and  PKP2 

by  irregularities  in  the  vicinity  of  the  Earth's  mantle-core 
boundary. 


1.  Introduction 


The  writer  (Haddon,  1972)  recently  drew  attention 
to  certain  inconsistencies  between  observational  data  on 
so-called  precursors  to  the  seismic  core  phase  PKIKP  and 
corresponding  theoretical  results  entailed  by  previous 
interpretations  of  the  precursors  involving  one  or  more 
transition  layers  surrounding  the  Earth's  inner  core.  He 
suggested  as  an  alternative  interpretation  that  the  pre- 
cursors may  originate  by  scattering  from  the  main  core 
phases  PKP^  and  PKP2  due  to  irregularities  in  the 
neighbourhood  of  the  mantle-core  boundary.  Subsequently, 
Cleary  and  Haddon  (1972)  examined  the  new  interpretation 
in  some  detail  and  assembled  a body  of  evidence  in  support. 
Further  support  has  recently  been  provided  by  Doornbos  and 
Husebye  (1972)  from  their  analysis  of  precursor  wavetrains 
recorded  at  the  Norsar  seismic  array  (see  Haddon  and  Cleary, 
1973) . More  recently.  King  (1973)  has  added  further  weight 
to  the  scattering  interpretation  with  results  from  his 
analysis  of  precursor  wavetrains  recorded  at  the  Warramunga 
seismic  array. 

In  addition  to  accounting  for  precursors  to  PKIKP, 
seismic  scattering  may  also  account  for  several  other  observed 
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seismic  "phases".  For  example,  Cleary  and  Haddon  (1973) 
have  suggested  that  much  of  the  body  wave  coda  following  P, 
including  the  so-called  precursors  to  PP  observed  by  Bolt 
and  colleagues  (1968)  and  others,  may  result  from  scattering 
in  the  crust  and  uppermost  part  of  the  upper  mantle. 

The  main  objection  to  the  scattering  interpretation 
raised  so  far  has  been  that  it  has  not  been  shown  that  the 
proposed  scattering  mechanism  can  account  for  the  observed 
amplitudes  of  the  phases  in  question.  It  is  therefore  of 
immediate  importance  to  investigate  the  mechanism  quantita- 
tively. The  present  paper  is  a first  step  towards  this  end. 

In  this  paper  a simple  theory  of  scattering  of 
elastic  body  waves  by  small  random  fluctuations  in  density 
and  elastic  parameters  is  developed  by  appropriately 
generalising  and  extending  certain  aspects  of  the  acoustic 
scattering  theory  given  by  Chernov  (1960)  . The  theory  below 
not  only  establishes  the  plausibility  of  the  proposed 
scattering  mechanism,  but  also  provides  an  adequate  model 
for  interpreting  observational  data  quantitatively.  For 
example,  application  of  the  theory  has  already  shown  that 
observed  amplitudes  of  precursors  to  PKIKP  can  be  accounted 
for  by  postulating  random  fluctuations  in  density  and  elastic 

parameters  of  order  one  per  cent  in  the  lowest  200  km  of 
the  mantle.  Further  numerical  details  relating  to-  the 
application  of  the  theory  will  be  published  in  a separate 


paper. 
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2.  Basic  Equations  for  Scattering 


Let  (i  = 1,2,3)  denote  the  rectangular 

cartesian  components  of  displacement  of  a point  P(x^)  in 
a perfectly  elastic  medium  in  which  the  density  p,  incom- 
pressibility k,  and  rigidity  y vary  slightly  from  their 
mean  values  from  point  to  point.  The  equations  of  motion 
for  small  displacements  in  such  a medium  may  be  written 


slr((k  - t “>6) 


+ 


+ 


3uk 

where  0 = -r — denotes  the  dilation. 
9xk 


(1) 


The  fluctuations  in  density  and  elastic  parameters 
will  be  denoted  by  Ap,  Ak  and  Ay,  so  that  p = pg  + Ap, 
k = kp  + Ak  and  y = Pg  + Ay,  where  pg,  kg  and  Pg  denote 
the  mean  values  and  it  is  assumed  that  | A p | <<  pQ,  | Ak | <<  kQ 
and  | Ay | <<  Pg.  For  the  present  pQ,  kQ  and  yQ  will  be 
assumed  to  be  constants.  (In  later  sections,  theory  developed 
on  this  basis  will  be  extended  to  cases  where  pQ,  kg  and  Pg 
are  slowly  varying  functions  of  distance  from  the  Earth's 
centre.)  Upon  substituting  for  p,  k and  y in  equation  (1) , 


we  obtain 
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32u. 


i /u  j.  1 x 36 

~ k f)  3 ^*0 

at"  0 3 0 


'0  - . 2 '"o  T I Mo'  ax.  ■ y0V  Ui 


a2u, 
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A • 1 

= - Ap  — 2^  + (Ak  + I AW>  ■§^7  + Au7"ui 

3 1 i 


(2) 


+ -g^-  (Ak  - y Ay)  6 + (Ay) 


3u.  3u. 

^7  + ^I 


When  Ap,  Ak  and  Ay  are  all  zero  the  equations 
(1)  and  (2)  reduce  to  the  usual  equations  of  motion  for  a 
homogeneous  medium.  Let  u^8  denote  any  solution  of 
equation  (2)  when  Ap,  Ak  and  Ay  are  all  zero  and  let 
ui  = u^  + ui1  denote  a corresponding  perturbation  solution 
of  equation  (2)  when  Ap,  Ak  and  Ay  are  non-zero.  Upon 
substituting  into  equation  (2),  assuming  that  1 ui 1 1 <<  | u^ 0 | , 
and  ignoring  terms  of  second  and  higher  orders  in  small 
quantities,  we  obtain 


a2u  1 


P°  " (k°  + ^ V Hr  - yov  V - 


Q^rft)  , 


(3) 


where 


el  _ 3uk 
8 " 7x7~ 
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r = (Xj) 
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and 


Ql  (r,t) 


~ 2 0 
3 u . 


- Ap  — —■  + (Ak  + ^ Ay)  96 
3t  J i 


,2.  0 


2 ' ' 1 L>UI  3x7  + ApV  ui 


+ W~  (Ak  ■ I Ap)0°  + 3!“  <A^> 
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Taking  the  divergence  and  curl  of  equation  (3)  gives 


and 


i - 7V 


“o  st 


r 1 

1 ~ - V2?1 

9 9 — V fc, 

V 9t 


= G(r,t)  , 


(5) 


= Hr,!)  , 


(6) 


where  r = (x^)  , 


p0a0  k0  + 3 U0  ' 


P0B0  u0  ' 


C1  = curl  (ui1)  , 0(r,t)  = div  (Qi)/(p0aQ2) 


= curl  (Qi)/(p0602)  . 


and 


T(r,t) 


The  resultants  of  the  secondary  scattered  waves 
originating  from  within  any  region  V of  the  medium  are 
given  by  the  following  solutions  of  equations  (5)  and  (6) 
(see,  e.g.,  Stratton,  pp.  424-428). 


0 


1 


0 (r' , t' )dV' 


(7) 


1 f 1 

Ttt  W 

. 

V 


'V  (r'  , t"  ) dV' 


(8) 


where  r'  = (x^1)  denotes  the  source  point,  r = (x^)  denotes 
the  field  point  at  which  the  solutions  are  to  be  evaluated 
at  time  t,  R'  = |r  - r' | is  the  distance  between  the 
source  and  field  points  and  t'  = t - R'/<*q  and  t"  = t - R'/Bq 
denote  retarded  times. 


The  equations  (5) , (6) , (7)  and  (8)  show  that 

under  the  influence  of  any  primary  wave  u^®,  each  element 
of  the  inhomogeneous  medium  becomes  effectively  a source  of 
scattered  waves  of  both  P and  S types.  In  the  following 
we  shall  restrict  attention  to  scattered  waves  of  P type 
originating  from  primary  waves  of  P type. 
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3.  P Wave  Scattering  from  a Primary  P Wave  by  Random 


Inhomogeneities 


In  general,  the  primary  P wave  disturbance  may 


be  represented  by 


9$ 


(i  ■=  1,2,3)  , 


where  4>  = 4>(r,t)  satisfies  the  wave  equation 


-i  2%  - v2*  - o . 
«o  3t 


Upon  substituting  from  equation  (9)  into  equation 


(4)  and  using  equation  (10)  , we  obtain 


Q.  = (Ak  + ^ " a02Ap)\iii 


+ (Ak  - 4 Ay)  .$  ..  + 2 (Ay)  .4>  . . , 
J f 1 » JJ  / J / 1 J 


where  the  subscripts  following  the  commas  denote  partial 
differentiations 


<t>  . . . 


9x^9x j ‘ 


We  shall  now  assume  that  the  fluctuations  Ap,  Ak 
and  Ap  are  given  by 

Ap/ PQ  = £H  , 

Ak/kQ  = mH  , 
and  Ap/pp  = nH  , 

where  l,  m and  n are  constants  and  H = H(r)  denotes  an 
isotropic  stationary  random  function  which  has  a correlation 
function  N given  by 


N = N (r)  = <H(r1)H(r2)>  , 


where  r here  denotes  | - r2 | . In  the  present  paper  we 
shall  further  restrict  consideration  to  the  particular  case 
where  N is  given  by 


N(r)  = exp  (-  r2/cr2)  , 


where  a is  a positive  constant  called  the  correlation 
distance  which  characterises  the  scale  of  the  random  fluctua 
tions  in  H.  It  may  be  noted  that  the  corresponding  root 
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mean  square  fluctuations  in  Ap,  Ak  and  Ay  are  given  by 


< (Ap;  > 


< (Ak) 


h 


£p0  ' 


mk0  , 


< (A y)  > = ny 


0 ' 


(15) 


so  that  H,  m and  n represent  the  magnitudes  of  the  root 
mean  square  proportional  changes  in  p,  k and  y. 

Upon  substituting  from  equation  (12)  into  equation 

2 

(11)  and  dividing  through  by  Pgao  • we  obtain 


-1 


Qi(p0a0  } 


= (Y1  - Y2)H*,ijj  + <Y1  - Y3)H,i*,jj  + Y3H,j*,ij  ' 


(16) 


where 


Yi  = (mk0  + £ np0)/(p0a02) 


m + ^-  (n  - m)  (3q/oiq)2  , 


^2  <lk0  + U ^p0^^p0a0  ) ~ ^ ' 


and 


Y3  2ny0/(pQa0  ) 


= 2n (0o/Qp)  . 


(17) 


Taking  the  divergence  of  equation  (16)  gives 


0(r,t) 


where 


Multiplying  corresponding  sides  of  equation  (7) 
by  their  complex  conjugates  and  averaging  over  the  ensemble 
of  possible  distributions  of  H gives  the  following  expres- 
sion for  the  mean  square  amplitudes  of  scattered  waves 
originating  inside  the  region  V. 


< I e 1 2>  = — 1 ~ <0(r'ft')0(r",t")>  dV'dV"  , (19) 

(4irr  R R 
V V 

where  r'  = (x^,)  and  r"  = (x^„ ) denote  the  integration 
variables,  R'  = |r  - r' | and  R"  = |r  - r" | denote  the 
distances  from  the  field  to  the  source  points,  t'  = t - R'/Oq 
and  t"  = t - R"/oiq  denote  the  retarded  times  and  the 
overbar  on  0(r",t")  denotes  the  complex  conjugate. 


Upon  substituting  from  equation  (18)  into  equation 
(19) , we  find  that  the  right-hand  side  of  the  resulting 
equation  can  be  expressed  as  a sum  of  integrals,  a typical 
member  of  which  is 


'23 


( 4 7T  ) 


a2a3 


FW7T 


R'R 


N 


t<J>,i'  j,  j' ] [ ^ , i ••  j«]dv'dV" 


(2 


V V 


where  N is  the  correlation  function  given  by  equation  (14) 
and  the  square  brackets  indicate  that  appropriate  retarded 
times  are  to  be  taken. 


Upon  transforming  the  integration  variables  in 
each  of  the  integrals  like  (20)  by 


2x j , = 2x j + ^ , 


2x j „ - 2x j - C j , 


j = 1,2,3 


we  obtain,  for  example, 


r r t 


23 


' l 


Go (x^,  £^,x^,t) 


’23 


R'  R" 


3 


(2 


1 


dx. dx2dx. 
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where,  by  equations  (14)  and  (21),  N is  now  given  by 


N = exp  (-  (Cx2  + C22  + Z>22)/°2]  ' 


(23) 


and  where  G23  ^xj,  • 'xi' denotes  the  transformed  form  of 

multiplied  by  a_a,.  In  equation  (22) 

fJ-JJ  # -*-  J ZJ 

the  domain  of  integration  V for  (x^,x2,x3)  is  identical 
with  the  domains  V for  (x^ , , x2 1 , x3 , ) and  (x^,, , x2„  ,x3„ ) 
while  the  domain  for  (£^,£^£3)  depends  on  (x^,x2,x3)  . 

However,  because  of  the  presence  of  either  the  factor  N or 
one  of  its  derivatives  in  the  integrands  of  each  of  the 
integrals  like  (20) , the  sum  of  which  comprise  the  right-hand 
side  of  equation  (19),  each  integrand  will  become  small  when 
|£|  becomes  large  compared  with  the  correlation  distance  o. 

It  follows  that  if  R' , R"  and  the  linear  dimensions  of  V 
are  all  large  compared  with  o,  then  a satisfactory  approxi- 
mation to  integrals  like  (20)  can  be  obtained  by  replacing 

the  factor  R'R"  by  R2,  where  R2  = { (x^  - x^)2  + (x2  - x2)2  + 
9 

(x^  - x^)  },  and  the  domain  by  the  infinite  domain  VB> 

We  thus  obtain,  for  example, 


23 


1 

ff 

( 4 7T  ) 2 j 

dx^dx2dx3 


9 3N 


2~r  G<xjl,£Jl,xi,t)d£1d£2d£3  , 


V 


( 4 TT  ) 


23  _ _ _ 

— Y dx^dx2dx3  , 

R 


say, 


(24) 
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I I 
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where  S. 


■ "V ■—  G(x,,Sp,x. 

3£.  35.  x 2 J 

i D 


Upon  combining  the  complete  set  of  terms  constituting 


the  right-hand  side  of  equation  (19),  we  obtain 


< er>  = 


( 4 it ) 2 R2  dxldx2dx3 


4 4 

where  S = 7 7 S 

L,  , mn 
m=l  n=l 


It  is  convenient  to  refer  to  S as  the  source  function  for 


scattered  waves. 


'?*■  * v • * •’ 

. s ' • 1 « 
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, A 

For  each  particular  point  r = r , say,  the 


functions  Smn  (above)  depend  (essentially)  only  on  the 
values  of  the  wave  function  $>  within  the  associated  region 


of  (Si'^2'^3)  where  N and  its  derivatives  are  significantly 
large.  We  shall  now  assume  that  within  each  such  limited 
region  4>  can  be  adequately  approximated  by  a plane  wave 
function  of  the  form 


4>  = 4>(n*(r  - r*)  + aQ(x  - t) ) , 


* 


! 

r, 


* * ★ 

where  r = (x^  , r = (xi  ),  n = n(r  ) denotes  a unit 

vector  in  the  direction  of  the  normal  to  the  wavefront  passing 

* * 

through  the  point  r , and  t = x (r  ) denotes  the  time  taken 

_* 

for  a wave  to  travel  from  its  origin  0 to  the  point  r . 


Now  let  x = Ty  denote  a transformation  of  coordi- 
nates from  the  x^  to  a new  "local"  rectangular  cartesian 

reference  frame  wh^-c^  has  its  origin  0^  at  the 

• , — ★ 

particular  point  x.^  = xi  in  x^-space,  its  y^  axis  in 
the  direction  of  the  vector  n and  which  contains  the  field 
point  P in  the  plane  y3  = 0.  Since  the  equations  (1)  to 
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(19)  are  all  independent  of  the  particular  cartesian 
reference  frame  used,  these  equations  apply  equally  well 
when  referred  to  a "new"  x^  coordinate  system  whose  axes 
coincide  with  those  of  the  y^  system  just  introduced.  In 
the  new  coordinate  system  equation  (26)  becomes 


<J>  = 4>(x1  + aQ(T  - t))  , 


(27) 


where  t is  the  time  taken  for  the  primary  wave  to  travel 

....  , ★ 

from  its  point  of  origin  0 to  the  origin  0^(x^  ) of  the 

new  coordinate  system.  For  the  remainder  of  this  section 

the  variables  x^,  xit,  x^„,  etc.  will  all  be  taken  to  refer 

to  the  new  coordinate  system,  unless  stated  otherwise. 


Corresponding  to  equation  (27)  the  equation  (18) 

reduces  to 


0(r,t)  = alHF"  + a2  |S-  F1  ♦ a3 

J-  dX-, 


jF+  a4V  HF  , 


(28) 


where  F = f(x1  + ocQ  (x  - t ) ) = 4>"  ^ + aQ(T  - t) ) denotes 
the  dilation  0Q  of  the  primary  wave  and  the  primes  on  F 
and  4>  denote  differentiations  with  respect  to  x^  Upon 
substituting  into  equation  (19)  and  using  certain  obvious 
symmetry  properties,  we  obtain 


? *-  * A.  : w't~'  V « . 
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+ 2a1.4»2‘MPl*P2 
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33N 


TFl’F2 


^ x ^ ii 

+ 2a2a4  Wf  ,22')Fl'F2 
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. 9xlt 
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?V12V22NF1F2JdV'dV"  , 


(29) 


where  Fx  = Ffs^,  F2  = F(s2),  F ^ ' = dF(s1)/ds1,  etc., 


= x^,  + R'  + cXq(t  “ t)  , 


S2  = Xl"  + R"  + ao  < 1 “ • 


and  R*  and  R"  again  denote  the  distances  between  the 
source  and  field  points. 

Upon  transforming  the  variables  (x^,)  and  (x^„) 
in  equation  (29)  by  the  transformation  (21)  and  approximating 
R'R"  by  R and  by  V^,  as  before,  we  obtain 


<|or> 


r?  dxidx2dx3 

K 
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2 3 ^ g4,  g ^ 

where  V denotes  *■  + *•  + *■  and  s.  and  s_  ari 

3£x  ac2  3^3^  1 2 

now  given  by 


= xx  + + R'  + aQ(T  - t)  , 


s-  = x,  - JsCi  + R"  + art(x  - t)  . 
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Although  the  equations  (25)  and  (30)  are  identical 

I in  form,  their  derivations  show  that  they  have  distinctly 

different  meanings.  For,  apart  from  referring  to  different 

coordinate  systems,  the  function  4>  given  by  equation  (26) 

has  been  assumed  to  approximate  the  general  function  <J> 

involved  in  equation  (25)  only  in  the  vicinity  of  the  partic- 

★ 

ular  point  x^  = x^  in  the  original  coordinate  system.  In 

respect  qf  the  functions  S given  by  equations  (25)  and  (31), 

however,  it  is  evident  that  under  the  assumed  circumstances 

the  value  of  the  function  S given  by  equation  (31)  at  the 

point  x^  = 0 in  the  new  coordinate  system  will  approximate 

the  value  of  the  function  S in  equation  (25)  at  the  point 
_ * 

x.^  = x^  in  the  original  coordinate  system  to  the  extent 
that  the  function  4>  given  by  equation  (26)  approximates 
the  general  function  4>  involved  in  equation  (25)  inside  the 
region  where  N and  its  derivatives  are  significant.  The 
equation  (31)  for  S evaluated  at  the  point  x-  = 0 there- 
fore provides  an  approximation  to  the  function  S in  (25) 

_ * 

at  the  point  x^  in  the  original  coordinate  system.  In  the 
next  section  we  shall  reduce  the  expression  (31)  for  S to 
a simpler  analytical  form. 


J 
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6 . Transformation  of  the  Source  Function  Integral 

Let  r'  denote  the  position  vector  O^P'  of  the 
source  point  P'(xit)  referred  to  the  local  coordinate 
system  introduced  in  the  previous  section,  R the  distance 
01P  from  the  origin  to  the  field  point  Ptx^,  n1  a 

unit  vector  in  the  direction  0][P  and  $ the  angle  between 
n and  n^,  as  shown  in  figure  1. 

Since  | r 1 | < < R and  O^P  is  in  the  plane  x^  = 0 
(by  choice  of  coordinate  system)  we  find  that  to  sufficient 
accuracy 


R’  = R - n1-r' 
and  similarly 

R"  = R - r^-r” 


R - x^,cos  <)>  - X2isin  4>  , 


R - x^„cos  <J>  - X2.(sin  <p  . 


(33) 


Transforming  these  expressions  by  use  of  the 
equations  (21) , then  substituting  into  equation  (32)  and 
putting  = X2  = x3  = 0 we  find  that 


si  = p?x  + q^2  + a0T ' * 

S2  = " P^i  " ^2  + a0T'  ' 


(34) 


and 
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2 

where  p = sin  <p/2  , 


q = - sin  <p/2  cos  <t>/2 


and  “q1'  = R + cIq(t  - t)  . 


Upon  integrating  with  respect  to  all  those 

terms  in  the  integrand  of  equation  (31)  involving  derivatives 
with  respect  to  we  find  that  all  such  terms  vanish 

because  and  F2  are  independent  of  £3  while  N and 

its  derivatives  vanish  at  the  limits  of  integration.  The 
equation  (31)  therefore  becomes 
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23 


where  = F(s^),  etc.,  and  and  s2  are  given  by 

equations  (34). 

Integrating  the  various  terms  in  equation  (35)  by 
parts  with  respect  to  and  £2,  as  appropriate,  and 

using  the  results  that 


q sfj"  (F1.F2)  " P ^7  (F1F2)' 


q2  <FiF2}  = P2  (FiF2)# 
H1  H2 


etc.,  we  obtain 


S = Re 


II  N{ai2W  + 2aia2  ifr  <W>  + 2ar 

V 1 


5 (Fi"F2) 
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2 3 
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‘2  2 '‘1  *2 
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+ a52  77T  <F1F2»  «1«2«3  • 

3C1 


(36) 


2 2 

where  a5  = (a3  + a^)  + a4q  /p  and  the  other  symbols  are  as 
previously  defined. 

After  some  further  manipulations  using  results  like 


Re 


N arq  (p1,p2")d51de2de3  - 


- Re 


N Wl  (Fl"F2,)d^1^2^3 


7 . Summary 


When  the  general  primary  wave  disturbance  4> 
(equation  (9))  can  be  adequately  approximated  by  plane  waves 
of  the  form  of  equation  (26)  within  each  subregion  of  V 
whose  linear  dimensions  are  of  order  a,  then  the  mean 
square  amplitude  of  scattered  waves  originating  inside  V 
and  arriving  at  the  field  point  P at  time  t is  given 
by  equation  (25)  where  the  source  function  S is  given  by 
equation  (37) . 


8 . Scattering  from  a Simple  Harmonic  Primary  Wave 


In  this  section  we  obtain  an  explicit  expression 
for  the  source  function  S (equation  (37))  for  the  case 
when  at  each  fixed  point  r = r of  the  scattering  region 
V the  primary  wave  can  be  represented  by 

Gq  = Aq  exp(iks)  , 

where  s = n*  (r  - r)  + aQ(x  - t)  and  where  the  unit  normal 
vector  n,  the  wavenumber  k and  the  amplitude  AQ  are 
assumed  to  vary  relatively  slowly  with  r,  and  t = t (r) 
denotes  the  time  taken  for  the  wave  to  travel  from  its 
source  to  the  point  P (r) , as  before.  For  this  case  we  have 

F(s1)  = Aq  exp(iks1)  , 

» 

F(s2)  = Aq  exp(iks2)  , 

where  s^  and  s2  are  given  by  the  equations  (34). 

Upon  substituting  from  equations  (23)  and  (40)  into 
equation  (37) , we  obtain 


When  |aq|,  r/  4>>  etc.,  vary  significantly  throughout 


the  region  V,  the  equation  (42)  would  generally  need  to  be 
evaluated  numerically.  If,  however,  equation  (39)  adequately 
represents  the  primary  wave  throughout  the  whole  of  V with 
k,  Aq,  and  n all  practically  constant,  and  if  also  R is 
large  compared  with  the  linear  dimensions  of  V,  then  the 
integrand  of  (42)  will  be  practically  constant  and  we 
immediately  obtain 

vr2 |a  | 2k4o3 

<|0|  > = 5 exp(-  k o sin  <)>/2)  , (43) 

16/irR^ 

where  V here  denotes  the  volume  of  the  scattering  region. 

For  the  particular  case  of  a fluid  medium  the 
equation  (43)  reduces  to  agreement  with  a corresponding 
result  given  by  Chernov  (p.  52) . 


The  observed  wavetrain  following  P frequently 
appears  to  have  a random  phase  character.  In  the  present 
section  we  therefore  consider  the  scattering  from  a primary 
wave  which  at  each  fixed  point  r = r of  the  scattering 
region  V can  be  represented  by 


J 

I I 

V 

j 

! 

k 

h 


0O  = AQF(s)  , (44) 

where  s = n- (r  - r)  + <xq(t  - t),  as  before,  AQ  denotes  the 
root  mean  square  amplitude,  and  where  F here  denotes  a 
real  stationary  random  function  with  a correlation  function 
M given  by 


M = <F(S;l)F(s2)  > = exp  (-  (Sl  - s2)2/X2)  . (45) 

i 

K ! I 

t In  this  equation  X denotes  the  correlation  parameter 

characterising  the  range  and  distribution  of  wavelengths 
in  the  signal. 

, Taking  the  statistical  average  of  equation  (37) 


over  the  ensemble  of  random  functions  F and  using  (45) , 
we  i n 


I 


- 31  - 


we  obtain 


<S>  = 


«V 


(3  - 48pu2/A2  + 64p2u4/X4) 


exp|-  u2  - (v2  + w2)/o2|dudvdw 


Integrating  (49)  gives 


<S>  = 


12r2A02(/?) 3a3 
,4fl  + sin2<t>/2 


and  substituting  into  equation  (25)  gives 


2 

<e“> 


r2_  2 3 

3 ( F A0  ° 

4/¥  V R2A4fl  + isi  , 


sin^4>/2 


As  in  the  previous  section,  if  the  integrand  of 
this  expression  varies  significantly  throughout  the  region 
V,  then  the  expression  would  generally  need  to  be  evaluated 
numerically.  If  the  integrand  remains  practically  constant 
throughout  V then  we  obtain 


<6  > = 


3vr2A02o3 

4 /ttR2 A 4 [ 1 + sin2 


sin2<)>/2 


where  V again  denotes  the  volume  of  the  scattering  region. 
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10 . Scattering  from  a Wavepacket  of  Finite  Length 

We  here  consider  scattering  from  a primary  wave 
disturbance  of  the  form 


0 0 «■  F (s)  = AqE  (s)  G (s)  , (53) 

where  s = n*  (r  - r)  + cxq(t  - t),  E(s)  represents  the 
envelope  shape  of  the  wavepacket  and  AQG(s)  represents 
either  the  simple  harmonic  wavetrain  given  by  equation  (39) 
or  the  random  wavetrain  given  by  equation  (44)  . Differen- 
tiating equation  (53)  gives 


F'  (s) 


1 + 


E'  (s)  7 (C  (s)1 

EjsT/  g (s)  J 


AqE(s)G'(s)  . 


(54) 


If  E(s)  is  a relatively  slowly  varying  function  compared 

with  G(s)  we  may  approximate  in  equations  such  as  (54)  by 

neglecting  the  derivatives  of  E(s).  For  example,  if 

E (s)  = e and  G(s)  = e1  s then  the  second  term  in 

the  brackets  in  equation  (54)  will  generally  be  less  than 
2 . 

, which  is  relatively  small  compared  with  unity  whenever 
the  wavepacket  contains  several  cycles. 

Under  the  assumed  circumstances  we  readily  obtain 
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< I 0 | 2 > = -i-  f 

16/?  J 


T2|AnE|2k4a3 


exp(-  k2c2sin2(j>/2)dV  , (55) 


<02> 


= -2-f 

A Sir  * 


r2|A0E|2a3 


4/?  ! R2 A 4 ( 1 + 


dV  , 


corresponding  to  the  equations  (42)  and  (51) , respectively, 
where  E = E(ctQT')  and  Og1'  = R + ag (T  " may  be 

noted  that  the  integrands  of  equations  (55)  and  (56)  are 
significant  only  in  those  parts  of  V where  E((XqT')  is 
significant.  The  scattering  region  may  therefore  be  taken 
infinite  when  appropriate. 


As  a final  example  of  the  application  of  equation 
(37) , we  consider  the  scattering  from  a primary  wave  pulse 
which  at  each  fixed  point  r = r of  the  scattering  region 
can  be  represented  by 


i i 


La 


.2  /,  2, 


e0  = Aq  exp  (-  s /A  ) , 


(57) 


where  s = n-  (r  - r)  + oiq(t  - t)  and  X here  characterises 


the  length  of  the  pulse.  In  this  case  we  must  evaluate 
equation  (37)  with 


2 / ■>  2 . 


F,  = F(s.)  = A_  exp  (-  s,/r)  , 


and  F _ = F ( s_ ) = An  exp(-  s„2/X2)  , 


(58) 


where  s.  and  s9  are  given  by  the  equations  (34). 


Consider,  for  example,  the  integration  of  the  first 
term  in  the  integrand  of  the  expression  (37). 





Upon  similarly  integrating  the  other  terms  in 


equation  (37)  and  combining,  we  obtain 


S = 


4 (/if) 

a4y5'2 


exp (-  2X2)  (g^4 


g2X 


g3} 


where  g^ 


4aiV 


'l 


g2  = ("  6ai2y2  + 2y(ci  - 3c3))  , > ( 

g3  = | (l  - y)2r2  + | c2y2  + (Cl  - | c2  + 3c3)y  , 


and  X and  Y have  the  same  meanings  as  in  equations  (60)  and 
(61)  . 


By  equations  (25)  and  (62)  the  resultant  mean  square 
amplitude  of  scattered  waves  arriving  at  P from  V is  given 
by 


<0 
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12 . Application  of  Scattering  Theory  to  Spherically 
Symmetrical  Earth  Models 


i 

i * 

r. 


It 


Consider  an  Earth  model  which  is  spherically 
symmetrical  except  for  small  random  fluctuations  in  the 
density  and  elastic  parameters  inside  a certain  region  V. 
Let  p0,  kQ,  Pq,  a0  and  denote  the  mean  density, 

incompressibility,  rigidity  and  P and  S velocities  at 
distance  r from  the  centre  of  the  model.  In  the  present 
section  we  shall  derive  a formula  for  the  resulting  P type 
scattering  when  a primary  P wave  propagates  through  the 
irregular  region  V.  We  shall  assume  that  the  primary  wave 
originates  from  a spherically  symmetrical  point  source  A 
and  that  ordinary  ray  theory  (see,  e.g.,  Bullen  (1963), 
chapters  7 and  8)  may  be  applied  to  calculate  travel  times 
and  amplitudes  of  both  the  primary  wave  and  the  secondary 
scattered  waves  excited  by  the  primary  wave. 


Let  I denote  the  energy  in  the  primary  wave 
emitted  per  unit  solid  angle  from  the  source  A and  let  eA 
be  the  angle  which  any  ray  makes  with  the  level  surface  r = r. 

ft 

through  A (see  figure  2) . Then  the  energy  E per  unit 
area  of  the  portion  of  the  wavefront  emerging  at  an  angle  en 
to  the  level  surface  r = rB  through  any  point  B is  given 
by  (Bullen,  p.  126) 


I 


Fig.  2 


The  primary  wave  ray  A B and  the  scattered  wave  ray  BC  for 
a primary  source  A,  a scattering  point  B and  a receiver  C in 
a spherically  symmetrical  earth  model:  x is  the  angle  of  inter 
section  of  the  diametral  planes  containing  AB  and  BC',  e , e , 
and  e are  the  angles  at  which  the  rays  intersect  the  level 
surfaces  through  A,  B and  C. 


where  AAB  denotes  the  angle  AOB  subtended  by  A and  B 
at  the  centre  O of  the  model  and  the  range  eA  to  eA  + deA 
corresponds  to  the  range  AAB  to  AAfi  + dAAfi. 


For  a simple  harmonic  wavetrain  the  amplitude  A 
of  the  wave  is  related  to  the  energy  E per  unit  area  of 
wavefront  by  (Bullen,  p.  128) 


where  the  factor  f depends  on  the  wavelength  and  period  of 
the  wave  and  the  density  of  the  medium.  A similar  relation- 
ship holds  generally  for  plane  waves  of  similar  waveforms. 

From  equations  (65)  and  (66)  we  obtain 
A 2 = f 1 C°S  eA  deA 

B B TT  sin  eB  sin  AAB  dAAB  ' 


where  An  is  the  amplitude  of  the  primary  wave  at  B and 

Jd  i 

fn  is  the  value  of  the  factor  f at  B.  The  expression  (67) 
will  be  applied  shortly  to  give  the  amplitude  of  the  primary 
wave  in  scattering  formulas  such  as  (55) , (56)  and  (64) . We 
shall  first  consider  effects  of  the  model's  velocity  structure 
on  amplitudes  of  scattered  waves. 


The  equations  (55)  (56)  and  (64)  can  all  be  written 


! in  the  form 

2>  = — L. 

(4n) 

where  AV  represents  an  infinitesimal  volume  element  of  V. 
The  resultant  of  the  scattered  waves  arriving  at  the  field 
point  P may  thus  be  regarded  as  the  sum  of  contributions 
originating  in  the  volume  elements  AV  comprising  V.  The 
contribution  from  a single  element  AV  may  be  represented  by 


< 0 


v SAV 

2 <•  2 ' 
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<A02> 


SAV  ‘ 
(4ttR)  2 


The  energy  per  unit  area  of  wavefront  arriving  at 
the  field  point  P from  AV  is  given  by 


AE 


SAV 

fB(4TTR)2 


This  contribution  may  be  regarded  as  coming  from  a point 
source  B within  AV  which  emits  SAV/ (f^x  (4tt)  units  of 
energy  per  unit  solid  angle  in  the  direction  BP.  The 
equations  (55) , (56)  and  (57)  all  apply  only  to  the  case 
where  Pg,  kg  and  vig  are  all  constants  throughout  the 


(68) 


(69) 


(70) 


region  connecting  B and  P.  We  shall  now  extend  these 
results  to  the  case  of  a spherically  symmetrical  model  where 
pQ,  kg  and  ^g  vary  with  r. 


If  a point  source  of  intensity  SAV/  (f0x  (4ft)  2) 
were  located  at, a point  B in  a spherically  symmetrical 
model  then  the  energy  AEg,  per  unit  area  of  wavefront 
originating  at  B and  passing  through  a field  point  C 
would  be  given  by  (cf.  Bullen,  p.  126) 


AE, 


SAV 


cot  ,eB' 


_ 2 sin  en  sin  Anr, 
(4?r)  fBrc  nfi  C . BC 


d2T 


BC 


dA‘ 


BC 


(71) 


where 


'B 


is  the  angle  which  any  ray  leaving  B makes  with 


the  level  surface  r = rQ  through  B,  ec  is  the  angle  which 


this  ray  makes  with  the  level  surface  r = rc  through  C, 


lBC 


is  the  travel  time  of  a wave  travelling  from  B to  C 


along  the  ray,  ABC  is  the  angle  BOC  subtended  by  B and 


nB  = rB/aB*  Thus  in 


C at  the  centre  0 of  the  model  and 
order  to  allow  for  geometrical  focusing  effects  associated 
with  velocity  structure  in  the  Earth  model  we  must  replace 
the  equation  (70)  by  the  equation  (71) . The  corresponding 
contribution  to  the  mean  square  amplitude  of  the  waves 
arriving  at  C is  then  given  by 


*.  Ji 


i •-  4 


42 
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d2T 


BC 
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BC 


(72) 


where  fc  is  the  value  of  the  factor  f at  C.  The 
resultant  mean  square  amplitude  at  C from  all  the  elementary 
sources  AV  then  becomes 


<0  > = 


(4 


it)2  /.  flr~7n 


fC  COt  sb' 


B‘C  "B  Sin  6C  Sin  ABC 


d2T 


BC 


dA" 


BC 


SdV  , 


(73) 


where  S denotes  the  source  function  appropriate  to  the 
particular  primary  wave  assumed  (e.g.  equations  (40),  (50) 
or  (62)). 

The  equation  (73)  will  now  be  referred  to  a coordi- 
nate system  which  is  particularly  suitable  for  numerical 
evaluation  of  the  integral.  Let  ip0  denote  the  angle 

D 

between  any  fixed  diametral  reference  plane  passing  through 
A,  and  the  diametral  plane  passing  through  the  focal  point 
A and  the  scattering  point  B.  We  shall  take  as  independent 
variables  \ p^,  eA  and  rfi.  (The  particular  advantage  of 
this  coordinate  system  is  that  the  surfaces  defined  by 

= ^1'  ^b  = ^2 ' eA  = el  and  eA  = e2  define  a tube  of 
rays  emerging  from  the  focus  A.)  In  the  coordinate  system 

taken  the  volume  element  dVB  bounded  by  the  surfaces  iJ>B, 

4>b  + eA'  eA  + deA'  rB  and  rB  + drB  is  given  bV 


..... --  ■ 


. fi  mm 
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13.  Scattering  from  Wavepackets  in  a Spherically  Symmetrical 
I Model 

For  scattering  from  a wavepacket  of  the  form  of 
equation  (53),  where  AqG(s)  is  given  by  either  equation 
(39)  or  (44),  the  equation  (76)  gives 


<|e|2>  = -Ip.1  f 

16 /nr,,  { 


C V, 


■ r2  [E (aQT ' ) ] 2k4a3  cos  e^  cot  e0’ 
J n0  sin  eR  sin  ec  sin  A0C 


d2T 


BC 


dA‘ 


BC 


exp(-  k2o2sin2<))/2)deAd^BdrB  , 


(77) 


and 


<02> 


JV 

4 /nr' 


r2[E(a0T,)l2o3 


V1  nB  A 


1 + 


4o 


~ sin  <)>/2j 


572 


cos  e^  cot  eQ' 
sin  efi  sin  e^sin  ABC 


2 

d^T 


BC 


dA 


BC 


de.dip-dr. 


(78) 


corresponding  to  equations  (55)  and  (56)  respectively,  where 
t'  is  now  given  by 


T ' = T + T - f- 

AB  iBC  ' 


(79) 


and  aQ  denotes  the  P velocity  at  the  point  (ea,^_,r  ). 

ABB 
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While  equations  (77)  and  (78)  are  in  an  ideal 
form  for  numerical  integration  in  some  applications,  in 
others  where  scattering  may  occur  at  or  near  the  lowest 
points  of  the  primary  wave  rays,  it  is  more  suitable  to 
change  from  the  variable  rg  to  sg,  where  sg  denotes 
the  length  of  the  ray  path  from  A to  B.  This  trans- 


formation may  be  accomplished  simply  by  substituting 


for 


sin  eB 


An  outline  of  a procedure  suitable  for  evaluating 


the  integrals  (77)  and  (78)  will  be  given  in  §15. 
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14 . Scattering  from  a Pulse-Like  Primary  Wave  in  a 
Spherically  Symmetrical  Model 

In  this  case  the  mean  square  amplitude  of  scattered 
waves  at  P is  given  by  substituting  from  equation  (62) 
into  equation  (76) . We  thus  obtain 


<e2> 


glx4  + 92x2  + g3  exP(“  2X2) deAd^BdrB  * 


(80) 


where 


J = 


V0' 


cos  eA  cot  eBf 


p-  r j - 

4/ttX^Y3/  r„  nn  sin  efi  sin  e,,  sin  A 


C "B 


BC 


2 

d T 


BC 


dA‘ 


BC 


(81) 


and  where  X is  now  given  by 


X = 


a 

T 


r(*. 


AB 


+ T 


BC 


- t) 


(82) 


In  general,  if  for  a particular  observation  point 
C and  a particular  time  t under  consideration,  the 
associated  values  of  the  scattering  angle  <t>  are  small  to 
moderate,  then  it  happens  that  X varies  only  relatively 
slowly  throughout  the  region  V.  In  such  circumstances  the 
integrand  of  (80)  would  therefore  vary  only  slowly  throughout 


V so  the  integral  can  be  evaluated  numerically  using  the 
same  procedure  as  for  the  integrals  (77)  and  (78)  (to  be 
outlined  in  §15). 


For  larger  values  of  4>,  however,  X will  change 

more  rapidly  and  the  main  contribution  to  the  integral  will 

. 2 

come  from  that  limited  part  of  V where  X <1.  Under 
these  circumstances  it  is  advantageous  to  integrate  over  X 
analytically  as  follows. 

We  first  change  the  integration  variables  in  (80) 
from  (eA'^B'rB)  to  {eA'^B'sB*  where  SB  denotes  the 
length  of  the  ray  from  A to  B.  The  expression  (QflL)  then 
becomes 

<02>  = ||  J gxX4  + g2X2  + g3] exp (-  2X2)sin  eBdeAd^fidsB  . 

Now  let  denote  a point  on  the  ray  defined  by  ^eA' ^B^ 

at  which  x'  = 0 (equation  (79)),  and  let  OjXyz  denote 

a local  rectangular  cartesian  coordinate  system  which  has 
its  origin  at  the  point  O^,  its  x-axis  in  the  direction  of 
the  ray  passing  through  0^  and  which  contains  the  field 
point  P in  the  plane  z = 0.  We  shall  suppose  that  the  ray 

through  O,  may  be  treated  as  coinciding  with  the  x-axis  in 


an  interval  -a  < x < b,  to  be  specified  shortly,  and  that 


in  the  same  interval  the  pencil  of  rays  defined  by  the 

deA  deA  di|»B  dip 

ranges  eA 5-  to  eft  + -3-  and  5-  to  ^ + -j- 

may  be  treated  as  parallel.  We  now  replace  dsn  by  dx  in 

equation  (82)  and  integrate  along  the  elementary  pencil  from 

x = -a  to  x = +b.  The  corresponding  contribution  to  the 


integral  (82)  is  then  given  by 


d<(T>  = 


x=b 

t f 4 2 

j[glX  + g2X2  + g. 

x=-a 


exp(-  2Xz)dx 


sin  eBdeAd^B 


Since  by  choice  of  0^  we  have  TAQ  + TQ  c - t = 0,  the 
equation  (82)  gives,  to  sufficient  accuracy. 


X = 


ao 

T 


(TAB  - taoJ  + tbc  “ To.c) 
= y|x  - (x  cos  <p  + y sin  <|>)j  , 


= 2 (px  + qy)/X  , 


2 

where  p = sin  <p/2  and  q = - sin  <J>/2  cos  <p/2  (cf.  equation 
(34)).  Changing  the  integration  variable  from  x to  X in 
equation  (83)  by  use  of  (84)  , we  obtain 


d<0  > = 


l^fg.X4  + g-X2  + g3]exp(-  2X2)dX 


sin  e^de-dif). 
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where  we  have  set  y = 0 in  equation  (84)  (since  the  pencil 
of  rays  is  indefinitely  thin) . Approximating  further  by 
treating  J sin  eB,  p,  g^,  g2  and  g3  as  constant  on  the 
interval  -a  < x < b and  integrating  (85)  with  respect  to 


X,  we  obtain 


d<02> 


K (a,b)  sin  e0  deAdipB  , 


where 


K (a,b)  = 


- 9,  — erf  + erf 


11  8x2^ 


23/2pa 


r23//2pb 


- a exp  (-  ^-)  + b exp  (- 


. l^(a3  expL  8rV]  + b3  expf 

,XJ  I X > 

/tt  £f23/2pa)  . jef23/2pb) 

+ g2j^37i[er£l— r^-J  + erfl-T^-| 


O 2.  2 

8p  b 


♦ erff^-ES 


' & a exp(-  ^T-)  + b exp(- 


l--*f23/2pa|  . _«f23/2pb 

+ g3  ^171 erf  l- r^-J  + erf[— IT" 


and  where  J sin  eB, 


at  the  point  0^. 


(87) 


p,  g^ , g2  and  g3  are  to  be  evaluated 
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So  far,  the  range  (-a,b)  has  not  been  defined. 

It  should  first  be  noted  that  no  matter  how  large  a and  b 
are  taken,  the  contributions  to  the  integral  (83)  from  outside 
the  range  -A/(2  sin24>/2)  < x < A/  ( 2 sin24i/2)  are  negligible. 

We  need  therefore  be  concerned  with  defining  a and  b only 
in  cases  where  it  is  appropriate  to  take  a or  b less  than 
A/ (2  sin  <}>/2)  . Such  cases  may  arise  when  the  region  V 
containing  the  irregularities  is  limited  in  extent  (as,  for 
example,  where  V is  a concentric  spherical  shell  r ^ < r < ^) . 
It  is  then  appropriate  to  take  a and  b so  that  x = -a 
and  x = b correspond  to  the  points  of  entry  and  exit  of  the 
ray  (eA, to  and  from  the  region  V.  Upon  combining  the 
contributions  (86)  from  all  the  ray  pencils,  equation  (80) 
becomes 


<Q2> 


’ 

I 


AJ  K(a,b)  sin  eR  de^di^ 
— 


(88) 


where  in  general  a = a(eA,i^B)  and  b = b(eA,i(jB)  and  E 
denotes  the  domain  of  (eA, corresponding  to  V^. 

If  the  dimensions  of  V and  the  magnitude  of  <J> 
are  such  that  2pa/A  and  2pb/A  exceed  unity,  then  we  may 
approximate  by  taking  a -*•  °°  and  b + °°.  Equation  (87) 
then  becomes 


T*;  C r 


1 
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T,  , _ v-  \ _ 3 /rr  _ ,1  /?  _ , /? 

K(a,b)  4 J572  gl  I ~pj2  g2  pT?  g3  ' 

which,  upon  substituting  from  equations  (36),  (38),  (41)  and 
(63) , becomes 


K (a,b)  = | T2  , (90) 

where  (remarkably!)  the  factor  T2  (equation  (41))  is  the 
same  factor  as  appears  in  the  scattering  formulas  for  simple 
harmonic  and  random  primary  waves. 


If  the  above  conditions  are  satisfied  for  all 
(eA,^B)  inside  any  domain  E(eA,^B) , then  by  equations  (81), 
(88)  and  (90)  the  scattering  from  that  part  of  V inter- 
cepted by  the  tube  of  rays  defined  by  E is  given  by 


<e2> 


3 

32/2  r 


fCI 


w 


r2a3 


cos  eA  cot  eB* 
E nBpA3Y5/2  sin  ec  sln  4bc 


d2T 


BC 


dA‘ 


BC 


deAdi(/B  , (91) 


where  the  point  B is  here  to  be  taken  as  the  point  on  the 
ray  (eA, ^3)  where  t1  = 0 (equation  (79)). 


■ fci*. . 


to  determine  TAfi,  AftB  and  efi. 

(5)  The  epicentral  distance  An„  from  the 
scattering  point  B to  the  observation  point  C is  then 
given  by 


cos  abc  = cos  A cos  aab  + sin  A sin  aab  cos  ^b  • 


(6) 

determine  the 
eB'  and  eQ 


Ordinary  ray  theory  is  then  applied  to 

travel  time  TB<-,,  the  angles  of  emergence 
2 2 

and  d TBC/dA  BC,  corresponding  to  ABC  and 


r 


(7)  The  angle  x between  the  diametral  planes 
containing  the  primary  wave  ray  AB  and  the  scattered  wave 
ray  BC  (see  figure  3)  is  given  by 


cos  x 


cos 


aab  cos 


sm  A 


AB 


Arc  - cos  A 
sin  ABC 


(8)  The  angle  <p  between  the  direction  of  the 
primary  wave  ray  at  B and  the  associated  scattered  wave  ray 
is  given  by 


cos  <p 


sin  eR  sin  e'  + cos  e cos  e ' cos  x • 
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(9)  The  integrands  of  equations  such  as  (77)  , (78) 
(80),  etc.,  are  thus  determined  numerically  at  each  point  of 
V and  hence  the  integrals  may  be  determined,  apart  from  the 
multiplicative  factor  (fcI). 

In  practical  applications  it  will  usually  be 
convenient  to  compare  the  amplitudes  of  scattered  waves 
with  the  amplitudes  of  the  associated  direct  waves  at  the 
Earth's  surface.  If,  in  equation  (67),  the  point  B is 
taken  to  coincide  with  a point  D at  the  surface,  then  we 
obtain 


<fDI>  cos  eA 

rD2  sin  eD  sin  AAD 


de. 


dA 


AD 


i 


(fpD  cot  eA 

^ sin  eD  sin  A 


AD 


AD 


AD 


Since  fD  = f^,  the  equations  for  scattered  wave  amplitudes, 
and  also  the  equation  (95)  for  direct  wave  amplitudes  all 
contain  the  factor  (f^I)  • The  relative  amplitudes  are 
therefore  independent  of  this  factor. 


A further  point  to  be  noted  is  that  if  either  the 
primary  wave  or  the  resulting  scattered  waves  pass  through 
interior  boundaries  at  which  discontinuities  in  pQ,  kQ  or 


I 


(95) 
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P0  occur,  then  the  amplitudes  should  be  modified  appro- 
/ priately  to  take  account  of  energy  losses  by  reflection  and 
conversion  from  P to  S.  In  some  applications  allowance 
for  attenuation  effects  would  also  need  to  be  made. 
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16 . Concluding  Remarks 

In  deriving  the  various  equations  given  in  this 
paper,  a number  of  approximations  and  simplifying  assump- 
tions have  been  made.  In  each  particular  application  it 
will  therefore  generally  be  necessary  to  consider  whether 
or  not  the  various  approximations  made  are  satisfactory. 

In  general,  it  is  to  be  expected  that  the  theory  will  give 
reliable  results  whenever  the  primary  wave  may  be  reasonably 
approximated  by  a plane  wave  in  each  part  of  the  scattering 
region  whose  linear  dimensions  are  comparable  with  the 
characteristic  size  of  the  postulated  inhomogeneities.  For 
short  period  waves  this  will  usually  be  the  case  except 
near  caustics  and  geometrical  shadow  boundaries.  Even  in 
these  exceptional  cases  it  is  not  unreasonable  to  expect 
on  physical  grounds  that  formulas  such  as  (77) , (78)  and  (91) 
will  provide  satisfactory  first  approximations  of  scattered 
wave  amplitudes.  In  this  connection  it  is  notable  that  in 
applications  involving  scattering  from  the  vicinity  of 
caustics,  the  integrands  of  equations  (77)  and  (78)  are  not 
singular  in  spite  of  the  infinite  amplitudes  associated  with 
caustics  on  simple  ray  theory.  In  fact,  equations  (77)  and 
(78)  show  that  on  the  theory  given,  the  scattering  amplitudes 
associated  with  any  particular  ray  tube  leaving  the  source 
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are  (almost)  independent  of  variations  in  the  cross-sectional 
area  of  the  ray  tube  along  its  length. 

The  above  theory  has  been  applied  to  calculate 
amplitudes  of  waves  scattered  from  PKP  waves  in  the  lowest 
200  km  of  the  mantle.  The  results  show  that  fluctuations 
in  density  and  elastic  parameters  of  order  one  per  cent  in 
that  region  fully  account  for  the  observed  amplitudes  of 
precursors  to  PKIKP.  Although,  as  indicated  above,  the 
results  on  scattering  from  the  vicinity  of  caustics  need  to 
be  treated  with  some  caution,  in  the  case  of  PKP  the 
numerical  calculations  give  practically  the  same  amplitudes 
for  waves  scattered  from  P before  entering  the  core  as 
for  waves  scattered  from  PKP  after  leaving  the  core.  The 
plane  wave  approximation  used  in  this  paper  should  be  fully 
satisfactory  in  the  former  case  at  least,  which  is  sufficient 
to  establish  the  plausibility  of  the  scattering  mechanism. 
Further  details  on  application  of  the  theory  to  the  PKIKP 
precursor  problem  will  be  published  in  a separate  paper. 
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